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VARIATIONS  IN  AIR  DENSITY  FROM  JANUARY  1972  TO  APRIL  1975 
AT  HEIGHTS  NEAR  200  km 

by 

Doreen  M.  C.  Walker 

SUMMARY 

Variations  in  air  density  have  been  determined  using  the  orbit  of  satellite 
Cosmos  462,  1971-106A,  which  entered  orbit  on  3 December  1971  with  an  initial 
perigee  near  230  km  and  inclination  65.75°,  and  decayed  on  4 April  1975.  Accurate 
orbits  determined  at  85  epochs  give  perigee  height  correct  to  about  200  m 
throughout  the  satellite's  lifetime.  Using  these  values  of  perigee  height  and 
orbital  decay  rates  from  NORAD  elements,  604  values  of  air  density  at  half  a scale 
height  above  perigee  have  been  evaluated.  These  densities  have  been  compared  with 
values  from  the  COSPAR  International  Reference  Atmosphere  1972,  taking  account  of 
variations  due  to  solar  activity  and  geomagnetic  disturbances,  and  day-to-night 
variations,  to  reveal  the  residual  variations  in  density  at  a series  of  standard 
heights,  245,  240,  232  and  213  km. 

The  main  residual  variation  is  semi-annual,  with  maxima  usually  in  April 
and  October,  and  minima  usually  in  January  and  July;  but  it  is  irregular  in  phase 
and  shape.  The  amplitude  of  the  semi-annual  variation  is  remarkably  constant 
from  year  to  year  between  1972  and  1975,  and  considerably  greater  than  that  given 
by  CIRA  1972'.  the  April/July  density  ratio  is  1.68,  not  1.32  as  in  CIRA',  the 
October-November  maxima  are  all  lower  than  the  April  maxima,  whereas  CIRA  gives 
the  opposite;  the  July  minima  are  18%  lower  than  the  January  minima,  as  opposed 
to  10%  in  CIRA. 

A standardized  semi-annual  density  variation  for  the  early  1970s  is  pre- 
sented, with  January  minimum  of  0.94,  April  maximum  of  1.28,  July  minimum  of 
0.77  and  October-November  maximum  of  1.22.  In  addition,  three  other  recurrent 
variations  are  recognizable:  in  each  year  the  density  has  a subsidiary  minimum 
in  May  and  maximum  in  June;  there  are  low  values  in  mid  November  and  high  values 
in  late  Decerober. 
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1 INTRODUCTION 

The  density  of  Che  upper  atmosphere  at  heights  between  ISO  and  1000  km  has 
been  studied  quite  thoroughly  from  the  analysis  of  satellite  orbits,  and  the 
results  of  these  studies  prior  to  1972  have  been  incorporated  in  Che  COSPAR 
International  Reference  Atmosphere  1972  . The  main  changes  in  density  are  due 
Co  the  effects  of  solar  extreme  ultraviolet  radiation,  the  day-to-night  varia- 
tion, geomagnetic  activity  and  the  semi-annual  variation.  Values  of  density 
obtained  from  orbital  decay  rates  of  satellites  can  be  cleared  of  Che  first 
three  effects  with  the  aid  of  CIRA  197 2 ^ by  utilizing  the  indices  of  solar 
activity  and  geomagnetic  disturbance  in  retrospect.  Removal  of  these  effects 
exposes  the  residual  variations,  a mainly  semi-annual  variation  with  maxima 

usually  in  April  and  October  and  minima  usually  in  January  and  July.  The  semi- 

2 

annual  variation,  first  detected  in  1961,  has  been  followed  until  early  1972  by 

3_8 

analysing  various  satellite  orbits  . The  aim  of  this  Report  is  to  evaluate 
the  density  at  heights  near  200  km  and  trace  the  semi-annual  variation  from 
January  1972  to  April  1975  using  the  satellite  Cosmos  462. 

Cosmos  462  was  launched  on  3 December  1971  and  was  one  of  a pair  of  satel- 

9 10 

lites  launched  by  the  Russians  to  test  high-speed  interception  * . Cosmos  462 

was  the  hunter  satellite  and  exploded  after  passing  the  target  satellite.  Cosmos 
459.  This  explosion  occurred  3.5  hours  after  launch  and  the  largest  piece  of 
the  satellite  remaining  in  orbit  was  designated  1971-106A. 

After  this  experiment,  I971-106A  had  an  orbital  period  of  about  105  minutes, 
perigee  height  230  km,  apogee  height  1800  km  and  inclination,  i,  65.7°;  Che 
satellite  remained  in  orbit  for  40  months,  without  further  disturbance,  and 
decayed  naturally  on  1975  April  4.90^^. 

During  its  lifetime  1971-106A  was  selected  for  high-priority  observing 

by  the  British  optical  and  radar  tracking  stations,  including  Che  Hewitt  cameras 

at  Malvern  and  Edinburgh.  Observations  were  also  available  from  the  United 

States  Navy,  kindly  supplied  by  the  US  Naval  Research  Laboratory;  from  the  Cape 

kinetheodolite  at  the  South  African  Astronomical  Observatory;  and  the  theodolite 

at  Jc'  Aen,  Finland.  Using  some  5400  observations,  accurate  orbits  at  85  epochs 

1 2 

have  mpuCed  at  RAE  throughout  the  satellite's  life. 

Tiie  ^/erigee  heights  obtained  at  the  85  epochs  have  been  used,  together 
with  the  orbital  decay  rates  calculated  from  NORAD  elements,  Co  derive  604  values 
of  air  density  at  heights  between  176  and  252  km. 
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2 ORBITAL  DECAY  RATE 

The  orbital  decay  rate,  n,  was  found  with  adequate  accuracy  from  the 
NORAD  elements  at  an  average  interval  of  two  days.  The  values  of  the  mean 
anomalistic  motion  n (rev/day)  were  differenced  to  give  values  of  An  , each  of 
which  was  divided  by  the  corresponding  time  interval  At  (day)  to  give  values 
of  n ■■  An/At  (rev/day^).  Any  erroneous  value  of  n causes  an  up-down  pattern 
in  n,  ie  one  value  too  high  the  next  too  low,  or  vice  versa.  After  inspecting 
the  values  of  n,  57  values  of  n responsible  for  up-down  patterns  were  elimin- 
ated; the  remaining  604  values  of  n are  plotted  as  circles  in  Figs  1 and  2. 

The  values  of  n from  the  85  RAE  orbits  are  plotted  as  triangles.  The  two  sets 
of  values  show  excellent  agreement,  especially  when  it  is  remembered  that  the 
RAE  values  are  for  specific  dates  and  the  values  from  NORAD  elements  are  averaged 
over  a time  interval  At. 

The  effects  of  solar  radiation  pressure  on  n are  negligible:  they  are 
calculated  to  be  less  than  0.02%  of  the  main  effects  due  to  air  drag. 

3 PERIGEE  HEIGHT 

3. 1 Variation  of  perigee  distance.  a(l  - e) 

Fig  3 shows  values  of  perigee  distance  a(I  - e) , plotted  as  circles,  from 
12 

the  85  orbits  determined  at  epochs  spread  quite  evenly  over  the  satellite's 
40-month  life,  where  a is  the  semi  major  axis  and  e is  the  eccentricity. 

For  most  satellites,  a(l  - e)  suffers  a regular  oscillation,  of  amplitude 
several  kilometres  and  period  a few  months,  caused  by  the  effects  of  the  odd 
zonal  harmonics  in  the  Earth's  gravitational  potential,  and  also  small  oscilla- 
tions due  to  lunisolar  perturbations.  Some  smoothing  of  the  observational  values 
of  a(I  - e)  is  required  and  in  the  past^  this  has  been  done  by  clearing  the 
odd-zonal  harmonic  oscillation  from  a(l  - e)  to  give  corrected  values,  which 
should  show  a steady  slow  decrease  as  a result  of  air  drag  (apart  from  small 
oscillations  due  to  lunisolar  perturbations).  A smooth  curve  was  then  drawn 
through  the  corrected  values  and  the  odd-zonal  harmonic  perturbation  restored 
to  give  smoothed  values  of  a(l  - e) . The  satellite  I97I-106A  has  an  orbital 

inclination  (65.7°)  at  which  the  oscillation  due  to  odd-zonal  harmonics  is  un- 
13 

usually  small  , its  amplitude  being  about  1.5  km,  and  the  period  is  unusually 
long,  about  two  years.  Consequently  the  values  of  a(l  - e)  show  only  a small 
and  slow  variation,  and  a smooth  curve  could  be  drawn  through  the  points  plotted 
as  circles.  However,  as  the  odd-zonal  harmonic  oscillation  is  so  small,  the 
oscillation  due  to  lunisolar  perturbations  is  of  the  same  order,  so  better 
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accuracy  can  still  be  achieved  by  calculating  the  zonal-harmonic  and  lunisolar 
perturbations,  removing  these  perturbations,  and  drawing  a smooth  curve  through 
the  corrected  values  of  a(l  - e) , Q say.  The  perturbations  due  to  zonal 


harmonics  and  lunisolar  effects  have  been  calculated  by  numerical  integration  at 

l4 

one-day  intervals  using  the  PROD  computer  program  . The  values  of  Q , which 


should  now  show  a steady  decrease  as  a result  of  air  drag  alone,  are  plotted  as 
crosses  in  Fig  3 and  a smooth  curve  has  been  drawn  through  the  points.  All  the 
values  have  standard  deviations  of  0.2  km  or  less. 


3.2  Perigee  height,  y 


The  actual  perigee  height,  y^,  above  the  Earth's  surface  is  found  by  first 


restoring  the  zonal-harmonic  and  lunisolar  perturbation  to  the  smoothed  values 


of  Q , giving  smoothed  values  of  perigee  distance,  a(l  - ®)smoothed’ 


then 


from  these  values  subtracting  the  local  Earth's  radius,  , at  latitude 
given  by 

R = R - 21.38  sin^(}> 

P P 


(1) 


where  R is  the  Earth's  equatorial  radius,  6378.14  km;  and  finally  adding  the 


small  amount  by  which  the  actual  path  of  the  satellite  departs  from  an  exact 
ellipse with  error  less  than  20  m,  we  may  write 


dr 


r 2 . 2 ' 

1 .388  1 — , ■ sin  u km 

1 + e 


(2) 


2 2 2 

where  co  is  the  argument  of  perigee.  Thus,  since  sin  (t>^  **  sin  i sin  u,  the 


values  of  perigee  height  are  given  by 

388  j^l  - pl^sin^coj 
- 6376.75  + ^17.76  - sin^u  . (3) 


2 2 

y - a(l  - e)  ..u  j ” 6378.14  + 21.38  sin  i sin  ui  + 1 . 
smoothed 


i(l  - e). 


smoothed 

The  values  of 

circles  in  Fig  4 and  a smooth  curve  has  been  drawn  through  the  points. 
4 EVALUATION  OF  AIR  DENSITY 


y^  calculated  at  the  85  epochs  using  equation  (3)  are  plotted  as 


4. 1 Determination  of  density 


The  air  density  p.  at  a height  iH*  above  perigee  was  obtained  from  n by 
16 

using  the  equation 


. 0.l57(1000«)  (J_-.  j I . ^ . 8.  . . 

ioV«  v»“*/  ttt;?  V 


2. 


c sin  i cos  2u 


(4) 


6 


c is  the  ellipticity  of  the  atmosphere  (taken  as  0.00335)  and  6 
16 


where 
area/mass  parameter 


IS  an 


17 


This  equation  was  originally  derived  with  the  aim  of  minimizing  errors 
in  due  to  errors  in  knowledge  of  scale  height:  H*  was  defined  as  the  'best 

estimate*  of  the  value  of  the  density  scale  height  H at  perigee,  and  the  equa- 
tion was  designed  to  allow  errors  of  up  to  25Z  in  H*  without  incurring  an  error 
of  more  than  1.2Z  in  p^.  Although  it  is  now  possible  to  assign  much  more 
accurate  values  to  H*  (correct  to  perhaps  5Z),  it  is  still  convenient  to  take 
advantage  of  the  insensitivity  of  p^  to  errors  in  H*  by  computing  density 
using  a constant  value  of  H*  over  time  intervals  during  which  H does  not 
vary  by  more  than  about  lOZ  as  a result  of  changes  in  perigee  height  and  solar 
activity.  For  1971-106A,  suitable  time  intervals  are  12,  12,  9 and  5^  months. 
During  the  first  of  these  chosen  time  intervals  the  variations  in  the  current 
value  of  H are  assessed  as  about  6Z  (sd)  due  to  variations  from  day  to  night 
and  variations  with  solar  activity,  and  about  4Z  (sd)  due  to  the  changes  in 
perigee  height  (±6  km) . Thus  the  total  error  in  H*  is  estimated  as  7Z  (sd) , 
giving  an  error  of  about  0.5Z  (sd)  in  p,.  The  other  three  intervals  have 

nearly  the  same  standard  deviations  in  p.. 

A 

The  constant  coefficient,  0.157,  outside  the  brackets  in  equation  (4)  is 

also  a function  of  the  gradient  of  scale  height  p (**  dH/dy),  and  the  numerical 

value  0.157  was  obtained^^  on  the  assumption  that  p ■ 0.1.  For  1971-106A  the 

appropriate  values  of  p are  (from  CIRA  1972)  between  0.13  and  0.18,  and  the 

departure  of  p from  its  asstuned  value  of  0.1  creates  an  error,  estimated as 

1.3Z  (sd) , in  the  coefficient  0.157.  However  this  error  in  the  formula  (4)  is 

automatically  cancelled  by  the  method  of  determining  6 (see  below)  and  need 

not  be  considered.  Further  errors  in  equation  (4)  are  caused  by  the  neglect  of 

^2  2 16 

small  terms  in  the  expansion,  which  are  of  order  e^,  ic^  cos  4a),  ^p^,  etc'  , 

2 

where  c ■ i)/H.  For  1971-106A  the  largest  of  these  terms  is 

2 

ic  cos  4(0,  which  has  a mean  value  of  0.010.  Thus  there  are  two  relevant  error 
terms  in  equation  (4),  0.5Z  from  variation  in  H,  and  1 .OZ  from  neglected  terms. 


leading  to  a total  error  in  p.  of  1.1 Z (sd)  for  given  values  of 


and  6 . 


The  values  of  H were  taken  from  CIRA  1972 ^ which  gives  the  variation  of 
H with  height  for  different  exospheric  temperatures.  Equation  (4)  is  valid  if 
3H*/a  < e < 0.2;  these  conditions  are  satisfied  until  three  weeks  before  decay. 

After  13  March  1975,  when  e < 3H*/a,  the  density  was  evaluated  using  the 
equation^ ^ 
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lOOOn* 


3iTa6n^  x 10® 


exp  (c 


cos 


2aj) 


1 + 2elj  (z*)/Iq  (z*)  + cI2(z*)/Iq(z*)  cos  2a)  ♦ c /4 


(5) 


where  z*  » ae/H*,  I (z*)  is  the  Bessel  function  of  first  kind  and  imaginary 

n 

argument  of  order  n and  is  given  by  Fig  7.4  of  Ref  16.  The  values  of 

3 

density  from  both  equation  (4)  and  equation  (5)  are  obtained  in  kg/m  if  H* 

and  a are  in  km,  the  area/mass  parameter  6 is  in  m /kg,  n is  in  rev/day 
2 

and  n in  rev/day  . 

The  values  of  air  density  from  equation  (4)  are  at  a height 

y.  ■ y + iH*,  and  the  values  p.  from  equation  (5)  are  at  a height 

A p A 

y^  “ y + XH*.  The  value  of  X is  0.5  when  ae/H*>  1 and  this  condition  is 
X p 

satisfied  for  all  the  values  of  p.  calculated.  So,  throughout  the  life  of  the 

A 

satellite  the  density  is  evaluated  at  a height 

The  value  of  the  area/mass  parameter  6 for  197I-106A  is  not  known  and  was 
chosen  to  give  values  of  density  in  conformity  with  CIRA  1972.  First,  the 
density  during  1972  was  calculated  from  equation  (4)  taking  a nominal  value  of 
0.01  for  6;  then  the  density  for  the  appropriate  local  time,  solar  activity, 
etc,  as  given  by  CIRAt  was  evaluated  at  six  dates  in  1972  when  the  semi-annual 
variation  was  near  its  mean  value  (February,  May,  September).  The  values  of  6 
needed  to  bring  the  values  of  density  from  equation  (4)  into  agreement  with 
those  from  CIRA  were  0.0069,  0.0074,  0.0072,  0.0080,  0.0069  and  0.0080  m /kg, 
giving  an  average  5 of  0.0074  with  rms  error  6%. 

4.2  Conversion  to  standard  height 

In  order  to  examine  day-to-day  variations  in  density  it  is  necessary  to 

convert  the  values  of  density  to  a fixed  height  y„  , say,  which  may  most 

c 

appropriately  be  taken  as  the  mean  value  of  y^  over  the  same  time  intervals 
as  those  used  in  section  4.1  for  H*  when  calculating  density  p^.  The  density 
pg  .at  height  y^  is  calculated  from  the  equation 


where  is  the  value  of  H at  height  y^.  Since  the  value  of  Hg  may  be  in 

error  by  about  7%  (sd)  due  to  day-to-day  variations  in  H,  and  the  rms  values  of 

(y^  " yg)  for  each  of  the  four  time  intervals  are  4.7,  5.3,  3.7  and  4.4  km 

respectively,  the  total  errors  in  p due  to  the  use  of  equation  (6)  will  be 

B 

1.4Z,  1.6Z,  1.4Z  and  1.6Z  respectively.  The  actual  time  intervals  taken  and  the 
values  of  the  various  parameters  are  given  in  Table  1 . 
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Table  1 

Four  time  intervals  and  values  of  parameters 


Time  interval 

Mean  y 

P 

H* 

yg  " ^A 

«B 

MJD 

km 

km 

km 

km 

41317  - 41682 
(1972) 

227 

36 

245 

39 

41683  - 42047 
(1973) 

223 

34 

240 

37 

42048  - 42320 

(1  Jan  - 30  Sep  1974) 

216 

32 

232 

35 

42321  - 42483 

(1  Oct  1974  - 12  Mar  1975) 

198 

30 

213 

32 

Between  the  end  of  the  fourth  interval  (13  March  1975}  and  decay  (4  April 
1975)  the  perigee  height  decreases  rapidly  and  conversion  of  to  a fixed 

height  is  no  longer  advantageous. 

5 PRESEHTATION  AND  ANALYSIS  OF  DENSITY 

5. 1 Variation  of  density  with  time 

The  values  of  density,  Pg,  obtained  at  fixed  heights  y^  for  each  of  the 
four  time  intervals  - 0245*  ^240*  ^232  *^213  ” plotted  as  circles 

joined  by  lines  at  the  top  of  Figs  5,  6,  7 and  8 respectively.  The  values  of 
Pj^  near  decay  are  shown  in  similar  format  in  Fig  9.  The  three  curves  below  the 
density  values  in  Figs  5 to  9,  in  descending  order,  are:  (1)  solar  radiation 
energy  on  10.7  cm  wavelength,  F^q  as  measured  by  MRC,  Ottawa;  (2)  sun-perigee 
angle,  with  the  local  time  at  perigee  in  hours  marked  above  the  curve  and  the 
latitude  of  perigee  in  degrees  marked  below  the  curve;  and  (3)  the  daily 
planetary  geomagnetic  index,  Ap,  as  given  by  the  Institut  fiir  Geophysik, 
Gottingen,  plotted  with  a 12-hour  time  lag. 

5.2  Comparison  with  CIRA  1972  model 

From  CIRA  1972,  values  of  density  can  be  obtained  for  the  same  times  as 
the  densities  evaluated  from  the  orbital  decay  rates.  The  CIRA  model  gives 
tables  of  density  against  height  for  a range  of  exospheric  temperatures.  To 
obtain  the  values  of  density  on  a particular  day,  the  exospheric  temperature 
has  to  be  evaluated  for  the  appropriate  conditions  of  (a)  solar  activity, 

(b)  local  time  at  perigee,  latitude  of  perigee  and  solar  declination,  and 

(c)  geomagnetic  disturbance.  By  comparing  c.  with  the  corresponding  value 
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from  CIRAt  the  residual  density  variations  (including  the  semi-annual  variation) 
can  be  traced.  Seasonal-latitudinal  variations  in  helium  have  no  effect  because 
the  proportion  of  helium  is  negligible  at  heights  near  250  km. 

5.2.1  Variation  with  solar  activity 

The  upper-atmosphere  temperature  and  density  reflect  changes  in  solar 
activity,  on  daily,  monthly  and  11-yearly  time  scales,  and  the  10.7cm  solar  flux, 
Fjq  is  generally  used  as  the  best  available  index  of  solar  EUV  radiation.  In 
order  to  evaluate  density  from  the  CIRA  model  at  a particular  height,  the  night- 
time minimum  temperature,  , in  kelvins,  is  computed  from  equation  (14)  in 
Part  3 of  CIRA 


Tc  ■ 1-94  F, 0.7* 


(7) 


where  ^]o  7 mean  flux  averaged  over  six  solar  rotations  (centred  on  the 

current  date)  and  7 bhe  flux  over  the  time  that  the  density  is  being 

^ determined  (usually  one  to  two  days).  The  value  of  T^  applies  for  an  undis- 

turbed geomagnetic  field  (Kp  ■>  0) . 

? 5.2.2  The  day-to-night  variation 

At  heights  above  about  200  km  the  density  of  the  upper  atmosphere  undergoes 
a diurnal  variation,  with  maximum  density  and  temperature  at  about  14  h local 
time  and  minimum  at  about  3 h local  time.  To  obtain  local  exospheric  temperature 
for  use  on  the  CIRA  model,  the  local  solar  time  at  perigee  and  latitude  of 
perigee  are  evaluated  at  intervals  of  approximately  20  days,  corresponding  to 
5°  changes  in  the  Sun's  declination;  then  the  appropriate  value  of  T^/T^  , the 
ratio  of  the  local  exospheric  temperature  at  perigee  T^  , to  the  minimum  night- 
time temperature  T^  , is  obtained  from  Table  1 in  Part  3 of  CIRA.  These  values 
of  are  plotted  against  time  and  a smooth  curve  is  drawn  through  the 

points.  Fig  10.  Values  of  can  then  be  read  off  the  curve  at  the  date  of 

each  density  value.  The  values  of  T^  thus  obtained,  using  T^  from  equation 
(7),  are  for  an  undisturbed  geomagnetic  field  (Kp  “0),  as  with  T^. 

5.2.3  Variations  with  geomagnetic  activity 

The  Sun  also  affects  the  upper  atmosphere  on  a day-to-day  basis  through  the 
incoming  particles  of  the  solar  wind.  When  the  Sun  erupts  with  flares,  intense 
streams  of  charged  particles  are  emitted  which  create  geomagnetic  storms.  These 
can  last  for  a few  hours  or  a few  days  and  cause  large  increases  in  density. 

The  geomagnetic  effect  on  the  temperature  in  kelvins  is  given  in  CIRA  1972  as 
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28Kp  + 0.03  exp(Kp) 


(8) 
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in  which  Kp  is  the  3-hour  geomagnetic  planetary  index.  In  this  study,  with 

density  evaluated  at  average  intervals  of  two  days,  the  daily  geomagnetic  index 

Ap  is  more  convenient  to  use  and  since  the  relationship  between  Ap  and  Kp 
18  . . 

is  known  , equation  (8)  gives  AT  as  a function  of  Ap  , and  is  shown  in  Fig  1 1 . 
As  Ap  ranges  from  2 to  400,  AT  varies  from  9 to  500  k'elvins.  Values  of  AT 
are  read  from  the  curve  of  Fig  1 1 for  values  of  Ap  corresponding  to  the  mean 
interval  over  which  the  density  is  being  estimated.  These  values  of  AT  are 
added  to  the  values  of  T to  give  the  appropriate  value  of  exospheric  tempera- 

ture,  T^,  corresponding  to  each  of  the  604  values  of  density  evaluated  from 
orbital  decay. 

5.2.4  Density  index,  D 

Using  the  values  of  exospheric  temperature,  T^,  corresponding  densities 
can  be  obtained  from  Table  6 in  Part  3 of  CIEA  1972.  These  values  of  density, 
^CIRA*  conditions  of  solar  activity,  diurnal  variation  and  geomagnetic 

disturbance  over  the  same  time  interval  as  the  densities  evaluated  from  the 
orbital  decay  rates. 

To  display  the  results,  the  values  of  density  from  the  orbital  decay,  p„, 

D 


are  divided  by  p. 


to  give  a 'density  index',  D , so  that 


for  each  of  the  four  fixed  heights,  y^,  given  in  Table  1.  The  values  of  D for 
each  of  the  four  time  intervals  are  plotted  as  circles  and  joined  by  straight 
lines  at  the  bottom  of  Figs  5 to  8. 

Between  13  March  1975  and  decay  (4  April  1975),  where  the  values  of 
density  are  not  converted  to  a standard  height,  the  values  of  D are  given 

by  ’^CIRA  evaluated  at  height  y^.  Thus  each  value  of  D 

is  for  a different  height,  but  still  represents  the  ratio  of  the  density  from 
orbital  decay  to  CIRA  density.  These  values  of  D are  plotted  at  the  bottom  of 
Fig  9.  The  values  of  y,  are  also  plotted  with  p,  at  the  top  of  Fig  9. 


RESULTS 


The  semi-annual  variation 


The  values  of  the  density  index,  D,  displayed  at  the  bottom  of  Figs  5 to 
9 give  a record  of  the  residual  variations  in  density  after  removal  of  the  main 
effects  of  solar  activity,  day-to-night  variations  in  density,  and  geomagnetic 
disturbances.  After  these  major  disturbances  due  to  the  Sun  have  been  accounted 
for,  the  most  important  remaining  variation  in  density  has  a semi-annual  recur- 


f*  T" 


rence  tendency,  generally  following  a pattern  of  maximum  density  in  April  and 
late  October,  and  minimum  in  January  and  July. 


In  the  results  from  1971-106A  the  course  of  this  semi-annual  variation  can 
be  traced  in  the  values  of  density  index,  D.  The  form  of  the  semi-annual 
variation  is  most  readily  revealed  by  drawing  two  parallel  curves  through  the 
jagged  observational  values,  as  in  Figs  5 to  8,  in  such  a way  that  the  majority 
of  the  observational  values  are  included  in  the  band  defined  by  the  two  curves, 
which  is  about  0.15  wide,  that  is  approximately  15%.  The  maximum  and  minimum 
values  in  the  semi-annual  variation  may  then  be  taken  as  the  maximum  and  minimum 
of  a curve  midway  through  the  band.  The  values  of  D thus  obtained,  with  the 
corresponding  dates,  are  given  in  Table  2. 

Table  2 

Semi-annual  maxima  and  minima  in  density  index  D 


Year 

Height 

January-February 

minimum 

March-April 

maximum 

July-August 

minimum 

October-Novembe 

maximum 

km 

Date 

1 

D 

I Date 

D 

Date 

D 

Date 

D 

1972 

245 

Feb 

20 

0.95 

Apr 

19 

1.30 

Aug 

2 

0.81 

Oct 

30 

1 .23 

1973 

240 

Jan 

18 

0.94 

Apr 

8 

1.29 

Jul 

17 

0.75 

Nov 

27 

1.18 

1974 

232 

Jan 

18 

0.91 

Mar 

18 

1.27 

Jul 

20 

0.74 

213 

Oct 

30 

1.26 

1975 

213 

Feb 

17 

0.97 

Mean 

values 

Feb 

3 

0.94 

Apr 

5 

1 .29 

Jul 

23 

0.77 

Nov 

8 

1.22 

CIRA 

values 

Jan 

16 

0.94 

Apr 

6 

1 .12 

Jul 

27 

0.85 

Oct 

27 

1 . 16 

Several  features  of  Table  2 and  Figs  5 to  8 call  for  comment: 

(1)  There  is  only  a small  variation  from  year  to  year  in  the  value  of  the 
density  index  at  a particular  maximum  or  minimum  of  the  semi-annual  variation: 
the  mean  values  are  within  5%  of  the  individual  values.  Three  of  these  mean 
values  however  differ  significantly  from  those  of  CIRA  1972;  the  April  and 
October  maxima  are  15%  and  5%  higher  than  in  CIRA  and  the  July  minimum  is  9% 
lower;  the  January  minimum  value  is  the  same  as  in  CIRA. 

(2)  The  October  maxima  are  all  lower  than  the  April  maxima,  the  mean 

value  being  5%  lower.  This  does  not  conform  with  CIRA,  where  the  October  value 

is  4%  higher  than  the  April  value.  The  CIRA  model  was  based  on  values  determined 

from  satellite  orbital  analysis  in  the  1960s,  when  the  October  maximum  was 
1 9 

usually  the  higher  . The  results  for  1972-75  do  however  agree  with  values 
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determined  from  1970-65D,  where  the  1970  November  maximum  was  4Z  lower  than  the 
1971  April  maximum  . The  high  April  maxima  may  therefore  be  regarded  as  a pheno- 
menon of  the  1970s. 

(3)  As  in  CIRA,  the  July  minimum  is  stronger  than  the  January  minimum, 
but  for  the  years  covered  by  this  study  the  July  minimum  is  18Z  stronger,  while 
CIRA  gives  it  as  only  10!S  stronger. 

(4)  The  shape  of  the  variation  is  not  sinusoidal,  and  the  maxima  and 
minima  are  not  at  regular  intervals.  The  most  striking  departure  from  a sinu- 
soidal curve  in  all  three  years  is  the  rapid  rise  to  a short-lived  maximum  in 
April,  followed  by  a long  slow  decline  to  the  flat  July-August  minimum.  The 
mean  dates  of  the  maxima  and  minima  are  given  in  Table  2 but  the  individual  dates 
differ  from  the  mean  by  as  much  as  19  days.  The  dates  of  the  CIRA  maxima  and 
minima  are  also  given  in  Table  2:  the  April  maximum  and  July  minimum  agree  to 
within  4 days,  but  the  January  minimum  and  October  maximum  differ  by  18  and  12 
days  respectively.  But  it  should  be  noted  that  although  the  mean  date  obtained 
for  the  April  maximum  agrees  with  the  CIRA  date  (April  5)  the  individual  dates 
differ:  the  maximum  in  1974  occurred  18  days  earlier,  on  March  18.  This 
emphasizes  how  variable  is  the  phase  of  the  semi-annual  variation. 

The  values  of  the  ratios  of  the  April  maximum  to  the  January  and  July 
minima,  and  the  October  maximum  to  the  July  and  following  January  minima,  are 
shown  in  Table  3.  The  values  for  each  individual  ratio  for  each  of  the  three 
years  are  suprisingly  constant  by  comparison  with  the  1960s  when  there  was  con- 
siderable variation  from  year  to  year.  The  factors  of  variation  are,  of  course, 
considerably  greater  than  those  from  CIRA^  shown  in  the  last  column  of  Table  3. 

Table  3 

Ratio  of  semi-annual  maximum/minimum  for  1972-1974 


1972 

72/73 

1973 

73/74 

1974 

74/75 

Mean 

CIRA 

April  maximum/ 
January  minimum 

1.37 

1.37 

1 .40 

1 .38 

1.19 

April  maximum/ 

July  minimum 

1.60 

1.72 

1.72 

1 .68 

1 .32 

October  maximum/ 
July  minimum 

1.52 

1.57 

1 ,70 

1 .60 

1.36 

October  maximum/ 
January  minimum 

1.31 

1 .30 

1.30 

1.30 

1.23 
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6.2  Other  variations  in  density  index.  D 

6.2.1  Day-to-day  variations 

The  jaggedness  of  the  density  index  in  Figs  5 to  8 reveals  how  the  density 
varies  from  day  to  day.  The  two  parallel  curves  drawn  to  show  the  course  of  the 
semi-annual  variation  enclose  most  of  the  jagged  observational  values.  So  the 
change  in  density  from  day  to  day  amounts  to  about  ±7%. 

There  is  very  little  likelihood  that  changes  in  cross-sectional  area  con- 
tribute to  the  jaggedness.  A satellite  that  rotates  many  times  while  in  the 
region  near  perigee  where  drag  is  important  (a  region  which  takes  about 
10  minutes  to  traverse)  will  have  an  effectively  constant  cross-sectional  area, 
the  average  over  a rotation.  The  flash  period  of  1971-106A  increased  from  about 
4 seconds  just  after  launch  to  about  20  seconds  early  in  1973,  with  no  abrupt 
changes;  during  1972,  therefore,  the  satellite  rotated  many  times  in  the  region 
near  perigee,  on  each  orbit.  So  the  jaggedness  in  1972  cannot  be  caused  by 
variations  in  cross-sectional  area,  and  since  the  appearance  of  the  graphs  for 
1973-75  is  very  similar,  there  is  no  reason  to  suspect  irregular  variations  in 
cross-sectional  area  during  1973-5. 

6.2.2  Effect  of  geomagnetic  storms 

The  observational  values  of  density  plotted  at  the  tops  of  Figs  5 to  8 
show  correlation  with  geomagnetic  activity.  The  density  index  D gives  the 
values  after  removal  of  these  effects  (see  section  5.2.3)  using  Jacchia's  equa- 
tion for  AT  in  CIRA  1972,  equation  (8).  However,  it  can  be  seen  that  on  some 
occasions  the  effect  has  not  been  fully  removed.  This  may  be  because  the 
geomagnetic  index  used,  Ap  , does  not  always  indicate  the  full  effect  of  a 
magnetic  storm  on  the  upper  atmosphere,  partly  because  of  the  incomplete  geo- 
graphical distribution  of  magnetic  observatories  whose  data  are  used  in  the 

20 

production  of  the  indices  , and  partly  because  these  measurements  of  effects 
at  ground  level  may  not  adequately  reflect  the  effects  on  the  upper  atmosphere. 
This  could  very  well  be  the  case  at  MJD  41438,  where  there  is  a sudden  increase 
in  density  but  only  a relatively  small  magnetic  storm  is  indicated  by  the  Ap 
inui.x  (Ap  • 40).  A similar  increase  occurred  at  MJD  42280  and  again  the  Ap 
index  was  small. 

It  may  be  that  the  geomagnetic  effect  could  be  more  adequately  removed  by 

representing  the  geomagnetic  variations  by  a latitude-dependent  approximation, 

. 2 1 
such  as  that  given  by  Roemer  , 
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AT  “ (21.4  sin  + 17.9)Kp  + 0.03  exp  Kp  kelvins  (10) 

where  ^ is  the  0.4-day  mean  of  the  3-hour ly  planetary  index  Kp  ; or  that 

22 

given  recently  by  Jacchia,  Slowey  and  von  Zahn 

AT  « A sin  <p  (11) 

where  A > 57.5Kp[l  + 0.027  exp  0.4Kp^  kelvins,  and  4>  is  the  latitude. 

Equation  (10)  gives  a larger  correction  than  equation  (8)  for  latitudes  greater 
than  30°,  and  would  therefore  be  more  satisfactory  for  some  of  the  values,  for 
example  at  MJD  41953  and  41985.  Equation  (11)  is  unacceptable,  however,  because 
it  gives  zero  effect  at  the  equator,  whereas  Fig  6 shows  that  the  storms  between 
41760  and  41780  are  accompanied  by  density  increases  of  about  20Z  and  30Z  res- 
pectively, although  perigee  is  very  near  the  equator. 

The  conclusion  drawn  is  that  the  Ap  index  does  not  always  adequately 
register  the  effects  of  magnetic  storms  on  the  upper  atmosphere.  However,  the 
formula  used  here,  equation  (8),  has  compensated  for  the  effects  reasonably  well, 
especially  so  in  the  major  storms  of  August  1972  and  July  1974. 

6.2.3  Further  recognizable  recurrent  variations 

Comparing  the  density  index  curves  for  the  years  1972-74,  after  allowance 
is  made  for  the  semi-annual  variation  and  any  residual  variations  at  times  of 
geomagnetic  disturbances,  other  features  which  recur  in  each  of  the  three  years 
become  apparent. 

Between  the  March-April  maximum  and  the  July-August  minimum  Figs  5 to  7 

reveal  a subsidiary  minimum  and  maximum,  with  the  minima  at  about  June  2,  May  25 

and  May  25  for  the  years  1972,  1973  and  1974  respectively,  and  the  maxima  about 

July  2,  June  15  and  June  10.  This  'June  revival'  in  density  also  seems  to  occur 

23 

in  1968-71:  from  analysis  of  1967-31A,  there  is  a minimum  near  May  25  and  a 
. . 4 

maximum  near  June  15  in  1968  and  a minimum  near  May  18  and  maximum  near  June  15 
in  1969;  analysis  of  1969-108A  indicates  a minimum  around  May  25  and  a maximum 
about  June  25  in  1970;  and  the  analysis  of  1970-65D,  indicates  a minimum  near 
May  29  and  maximum  near  June  28  in  1971.  So  the  'June  revival'  shows  itself  in 
seven  successive  years,  with  minimum  and  maximum  around  May  26  and  June  20  res- 
pectively. There  is  also  some  evidence  of  this  'June  revival'  in  earlier  years 

24 

from  the  data  in  1966-69  presented  by  Wulf-Mathies  and  in  the  average  values 

1 9 

for  the  1960s  of  Jacchia  and  Slowey 
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Two  more  similarities  are  worth  a mention:  first,  there  is  a group  of 
values  below  the  semi-annual  band  from  14-23  November;  second,  another  group  lies 
above  the  semi-annual  band  from  19-27  December.  This  feature  is  apparent  in  all 
three  years  and  also  in  the  results  from  1969-20B^  for  1969.  Results  for  the 
years  1971  and  1972  are  not  available  in  sufficient  detail  to  determine  if  the 
same  trend  is  present. 

6.3  Variation  of  density  with  height 

The  variation  of  density,  , with  height  is  given  in  Fig  12.  For  dates 
before  13  March  1975  (when  y > 208  km)  the  values  of  p^  plotted  are  those  near 
the  averages  of  both  the  semi-annual  variation  and  the  day-to-night  variation, 
defined  as  the  values  of  p^  when  0.95  < D < 1.05  and  1.12  < T^^/T^  < 1.18  , 
and  also  when  geomagnetic  disturbances  are  small,  Ap  < 30  . After  13  March  1975 
all  the  points  are  plotted  except  those  with  Ap  > 30  ; the  points  plotted  are 
values  of  p^  , because  e < 3H*/a  (see  section  4.1). 

The  values  of  density  are  plotted  with  different  symbols  to  indicate  into 
which  bands  of  exospheric  temperature  they  fall.  The  three  curves  drawn 
represent  the  CIRA  1972  values  of  density  for  exospheric  temperatures  of  800, 

900  and  1000  kelvins. 

Fig  12  shows  there  is  good  agreement  between  the  observational  points  and 
the  curves  of  CIRA  1972.  Agreement  at  heights  near  245  km  is  to  be  expected, 
because  the  value  of  6 was  chosen  to  give  densities  consistent  with  CIRA  at 
that  height.  Fig  12  therefore  shows  that  CIRA  1972  gives  the  correct  variation 
of  density  with  height  between  250  km  and  175  km. 

7 DISCUSSION 

7. 1 Solar  activity  and  diurnal  variations 

Inspection  of  Figs  5 to  8 reveals  that  the  variations  due  to  solar  activity 
and  diurnal  variations  have  been  successfully  removed. 

The  values  of  ^|q  7 have  been  used  (section  5.2.1)  to  remove  the  varia- 
tions due  to  solar  activity  from  the  p^  values,  and  the  graph  of  density  index 

O 

shows  no  correlation  with  ^jq  7 • For  example  in  Fig  5 at  MJD  41360  there  is 
a steep  rise  in  P245  ^ peak  at  MJD  41370,  corresponding  to  a similar  rise 

in  Fjq  ^ ; a look  at  the  values  of  D shows  that  this  increase  has  been 
successfully  removed. 

The  diurnal  variations  give  a maximum  in  density  around  14  h local  time 

and  a minimum  at  03  h local  time.  This  variation  is  visible  in  p,^  but  no 

J5 
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trace  of  it  (or  over-correction  of  it)  is  discernible  in  the  values  of  D , so 

it  can  be  assumed  that  this  variation  is  correctly  modelled  in  CIRA  1972.  The 

amplitude  of  the  exospheric  temperature  day-to-night  variation  is  taken  as  0.3 

25 

in  CIRA  1972,  but  Broglio  et  dl  have  suggested  that  this  value  should  vary 
with  altitude  and  give  a value  of  0.19  at  240  km  height  decreasing  to  0.14  at 
210  km  height;  however,  there  is  no  indication  that  this  change  is  needed  in 
the  present  study. 

7.2  Effects  of  geomagnetic  disturbance 

The  increases  in  density  at  times  of  geomagnetic  disturbance  have  been 

allowed  for  by  using  equation  (8)  (see  section  5.2.3).  The  increases  in  density 

have  not  always  been  completely  removed  and  various  reasons  are  possible. 

26 

Nisbet  et  dl  suggest  that  the  effects  of  geomagnetic  activity  are  stronger  in 

the  morning  than  the  afternoon.  At  the  times  of  six  geomagnetic  storms, 

HJD  41438,  41865,  41985,  42280,  42335  and  42360,  the  effects  of  geomagnetic 

activity  have  not  been  completely  removed  from  the  density  index  values,  but  as 

three  of  these  storms  occurred  in  the  morning  and  three  in  the  afternoon  the 

suggestion  of  Nisbet  et  dl  does  not  seem  to  apply  here.  A second  possibility 

is  that  equation  (8)  should  include  some  effect  of  latitude,  but,  as  already 

explained  in  section  6.2.2,  the  use  of  latitude-dependent  equations  is  unhelpful. 
27 

Trinks  et  dl  conclude  that  the  global  geomagnetic  index  is  inappropriate  for 
correlation  with  locally  confined  data.  This  confirms  the  conclusion  of 
section  6.2.2  as  the  most  likely  reason  for  the  incomplete  removal  of  the  storm 
effects,  as  all  the  values  of  density  calculated  here  are  for  a localized 
section  of  the  orbit  - an  arc  of  40°  centred  on  perigee. 

7.3  Year-to-year  variations  in  semi-annual  effect 

From  the  analysis  of  satellite  orbits  during  the  1960s  a three-year 

oscillation  in  the  strength  of  the  semi-annual  variation  was  detected. 

28 

King-Hele  and  Walker  first  suggested  that  the  amplitude  of  the  semi-annual 

variation  oscillated  with  a period  of  about  33  months,  with  maxima  in  early 

3 29 

1959,  late  1961,  and  1964;  and  Cook  * showed  that  a three-year  periodicity 

30 

existed  from  1958  to  1967.  Voiskovskii  et  al  give  semi-annual  variations  at 

240  km  height  derived  from  analysis  of  the  orbits  of  Cosmos  satellites,  and 

their  amplitudes  have  clear  maxima  in  1964  and  1967;  but  for  the  years  1968-70 

31 

the  picture  is  less  clear.  Marov  and  Alpherov  give  results  from  Cosmos 

satellites  from  1961  to  1969,  which  again  show  maxima  in  1964  and  1967. 

32 

Wulf-Mathies  et  dl  give  the  results  from  analysis  of  Explorers  19,  24  and  39, 
which  also  show  a large  semi-annual  effect  in  1967. 
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In  Fig  13  the  results  given  in  Table  3 are  combined  with  others  from  1966  ■ 

4—8  33—35 

to  1972  Mean  values  have  been  taken  of  the  ratio  of  the  October  maxi- 

mum to  the  following  January  minimum  and  April  maximum  to  the  same  January  mini- 
mum each  year,  and  the  values  are  plotted  as  triangles  joined  by  a dashed  line. 

Similarly  mean  values  of  April  maximum  to  July  minimum  and  October  maximum  to 
July  minimum  each  year  are  obtained  and  plotted  as  circles  joined  by  a dash-dot  i 

line.  The  full  line  in  Fig  13  is  the  mean  of  the  two  curves  already  mentioned  | 

and  should  represent  the  year-to-year  changes  in  the  amplitude  of  the  semi-annual  | 

variation.  The  1967  maximum  is  clearly  visible  but  there  is  no  maximum  in  1970  | 

as  would  be  expected,  and  after  1971  the  value  remains  almost  constant  with  | 

average  value  1.48  ±0.04.  > 


So  the  three-year  oscillation  apparent  in  the  1960s  disappeared  during  the 
first  half  of  the  1970s. 

7.4  Standardized  semi-annual  variation  for  1972-4 

The  results  show  a semi-annual  variation  which  differs  considerably  from 
that  of  CIEA  1972  and  also  agrees  with  the  results  in  1970-71  ^ that  the  April 
maximum  was  stronger  than  October  maximum.  So  it  is  worth  trying  to  construct  a 
standardized  semi-annual  variation  in  density  for  1972-74  at  heights  of  200- 
250  km. 

Fig  I4a  shows  the  variation  of  D , the  value  of  the  density  index  D 
given  by  curves  midway  between  the  dashed  curves  in  Figs  5 to  8.  The  choice  of 
a standard  curve  for  D is  difficult,  but  as  the  individual  values  at  each 
particular  maximum  and  minimum  are  within  5%  of  the  mean  values,  the  standard 
curve  has  been  constructed  by  adopting  the  mean  maximum  and  minimum  values  given 
in  Table  2 and  choosing  dates  which  are  acceptable  for  the  majority  of  the  years. 
Fig  14b  shows  this  standard  curve:  the  dates  above  the  curve  are  the  mean  dates 
of  the  maximum  and  minimtun  given  in  Table  2.  The  actual  dates  from  the  standard 
curve  for  the  January  minimum,  April  maximum,  July  minimum  and  October-November 
maximum,  given  below  the  curve  in  Fig  14b,  are  January  24,  April  5,  July  21  and 
November  5.  As  it  happens  the  only  major  departure  from  the  mean  dates  given  in 
Table  2 is  that  for  the  January  minimum;  this  is  ten  days  earlier  than  the  mean 
date  of  Table  2,  and  gives  more  weight  to  the  results  for  1973  and  1974  obtained 
here,  the  date  for  1971  ^ (January  5)  and  the  date  (January  16)  for  the  1960s 
given  by  CIRA  1972  * . Although  the  recommended  curve  does  not  apply  for  any 
single  calendar  year,  it  is  close  to  the  variation  during  the  twelve  months 
starting  on  1 September  1972. 
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Fig  14b  is  offered  as  a method  of  determining  density  in  the  early  1970s: 
first,  density  should  be  calculated  from  CIRA  1972 » taking  account  of  the 

appropriate  local  time,  latitude,  solar  activity  and  geomagnetic  disturbance 
but  ignoring  the  CIRA  semi-annual  variation.  Then  should  be  multiplied 

by  the  standardized  D of  Fig  14b.  The  results  apply  for  heights  of  200-250  km, 
but  the  semi-annual  effect  is  not  strongly  dependent  on  height  and  the  results 
may  be  valid  down  to  about  150  km  and  up  to  300  km.  As  the  actual  semi-annual 
variation  differs  from  that  of  CIRA,  there  is  no  reason  why  the  mean  value  of 
the  standardized  D over  a year  should  be  1.0:  the  actual  mean  value  in  Fig  14b 
is  1.030. 


7.5  Possible  explanations  for  the  semi-annual  variation 

Many  theories  for  the  origin  of  the  semi-annual  variation  have  been  pro- 

3 

pounded,  but  none  of  them  offers  a con^lete  explanation.  In  1969  Cook  discussed 
in  detail  numerous  explanations  proposed  in  the  1960s,  and  enqihasized  the  con- 
tinuation of  the  effect  down  to  about  90  km;  but  he  concluded  that  the  semi-annual 
effect  should  be  regarded  as  "one  of  the  least  understood  phenomena  of  the 
atmosphere" . 

The  name  'semi-annual  variation'  describes  an  effect  that  really  comprises 
two  variations  - one  annual  and  the  other  (and  larger)  semi-annual.  An  annual 
variation  arises  because  of  the  variation  in  the  Earth's  distance  from  the  Sun; 
the  solar  radiation  received  in  January  is  7%  greater  than  in  July,  so  the 
density  in  January  at  heights  near  240  km  is  expected  to  be  about  6Z  greater 
than  in  July.  This  accounts  for  much  of  the  difference  between  the  January  and 
July  values,  but  the  semi-annual  variation  remains  to  be  explained. 

Of  the  various  explanations  of  the  semi-annual  effect  that  have  been 

proposed,  the  most  promising  seem  to  be  those  which  attribute  the  effect  to 

seasonal-latitudinal  variations  in  the  mesosphere  extending  into  the  thermosphere. 
36 

Volland  compared  the  effects  produced  by  heating  at  various  heights  in  the 
thermosphere  and  concluded  that  the  variation  of  the  semi-annual  effect  with 
height  was  best  modelled  by  heat  input  at  the  base  of  the  thermosphere.  He 
suggested  that  a leakage  of  wave  energy  flux  of  a tidal  wave  into  the  thermo- 
sphere from  the  lower  atmosphere,  due  to  changes  in  mesospheric  wind  systems 

which  are  often  semi-annual  in  character,  produces  the  observed  density  variations. 
31 

Marov  and  Alpherov  conclude  that  the  semi-annual  density  variations  in  the 

thermosphere  are  linked  with  seasonal  circulation  of  the  lower  atmosphere.  Ching 
. 37 

and  Chiu  point  out  that  the  total  solar  radiation  impinging  on  the  Earth,  in 

38 

both  the  lower  and  upper  atmosphere,  has  a semi-annual  variability.  Fukuyama 
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finds  a semi-annual  variation  in  night  airglow  intensity,  which  supports  the 

39 

view  of  linkage  between  lower  and  upper  atmosphere.  Groves  has  described  semi- 
annual variations  in  density  in  the  100-km  region  which  link  with  the  thermo- 

40  . . 

spheric  semi-annual  variation.  Groves  has  also  investigated  the  absorption  of 
solar  radiation  by  water  vapour  and  ozone,  and  concluded  that  the  observed 
variations,  particularly  the  low  value  in  July,  agree  with  the  thermospheric 
semi-annual  variation. 

7.6  Congjarison  with  Jacchia's  new  models  (1977) 

41 

After  the  completion  of  this  work,  Jacchia's  latest  atmospheric  models 
were  issued.  He  states  that,  in  selecting  the  form  and  amplitude  of  the  semi- 
annual variation,  "we  still  prefer  to  use  the  model  of  J71",  which  is  the  same 
as  the  CIBA  1972  model.  However,  the  tables  he  gives  show  a slightly  increased 
aiiq>litude:  the  values  at  the  January  minimum,  April  maximum,  July  minimum  and 
October  maximum  for  comparison  with  those  of  Table  2 are  0.94,  1.14,  0.83  and 
1.18  respectively.  Jacchia's  new  values  are  therefore  generally  nearer  than 
those  of  CIRA  1972  to  the  variations  found  here  for  1972-5,  but  the  airq)litude 
is  still  too  small:  in  particular  his  ratio  of  the  April  maximum  to  the  July 
minimum  is  1.37,  as  opposed  to  1.68  found  here. 

8 CONCLUSIONS 

Analysis  of  the  orbit  of  1971-106A  during  its  40-month  lifetime  has  yielded 
604  values  of  air  density  at  heights  near  200  km. 

Variations  in  the  density  values  due  to  solar  activity  and  diurnal  varia- 
tions have  been  successfully  removed  using  CIRA  1972.  Allowance  was  also  made 
for  the  increases  in  density  at  times  of  geomagnetic  disturbances  using  CIRA  1972 
(see  section  5.2.3),  but  this  was  not  always  so  successful  and  various  possible 
reasons  have  been  given  (section  7.2).  After  removal  of  these  effects  and  con- 
version of  density  to  a standard  height,  a record  of  the  residual  variations  in 
density  between  January  1972  and  April  1975  emerges. 

The  main  residual  variation  has  a semi-annual  pattern  and  the  maximum  and 
minimum  values  of  t'he  variation  are  given  in  Table  2.  The  results  show  that: 

(1)  There  is  only  a small  variation  from  year  to  year  in  the  value  of  the 
density  index  at  a particular  maximum  or  minimum  of  the  semi-annual  variation; 
but  three  of  the  mean  values  are  significantly  different  from  those  of  CIRA  1972, 
always  giving  a stronger  variation  than  CIRA. 
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(2)  The  October  maxima  are  all  lower  than  the  April  maxima,  the  mean  value 
being  SZ  lower.  This  does  not  conform  with  CIRA  1972  where  the  October  maximum 
is  4Z  higher  than  the  April  value. 

(3)  As  in  CIRA  1972,  the  July  minima  are  stronger  than  the  January  minima,  but 
18Z  stronger  where  CIRA  gives  lOZ. 

(A)  The  shape  of  the  variation  is  not  sinusoidal  and  the  maxima  and  minima  are 
not  at  regular  intervals. 

(5)  When  the  ratio  of  the  April  maximum  to  January  and  July  minima  and  the 
ratio  of  the  October  maximum  to  the  July  and  following  January  minima  are  com- 
pared (see  Table  3) , the  values  for  each  individual  ratio  for  each  of  the  three 
years  are  surprisingly  constant,  in  contrast  to  the  1960s,  when  there  were  wide 
variations  from  year  to  year.  The  ratios  of  maximum/minimum  density  in  1972-5 
are  greater  than  those  of  CIRA  1972,  and  close  to  the  average  of  the  maximum 
values  of  the  ratios  in  the  1960s  (see  Fig  13). 

The  variation  of  density  with  height  between  176  and  252  km  obtained  from 
analysis  of  1971-106A  is  in  close  agreement  with  the  CIRA  1972  model  (see  Fig  12). 

In  addition  to  the  main  semi-annual  variation,  the  air  density  shows  three 
other  recurrent  features.  The  most  in^ortant  is  a minimum  in  late  May  and  a 
maximum  in  mid  June.  There  is  also  a mid-November  minimum  and  a late-December 
maximum  (see  section  6.2.3). 

Since  the  amplitude  of  the  semi-annual  variation  in  1972-5  is  so  consistent, 
a standardized  curve  for  the  semi-annual  variation  in  the  early  1970s  is  pre- 
sented (Fig  14),  applying  for  heights  of  200-250  km. 
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Fig  1 Orbital  decay  rate  n for  Cosmos  462,  from  January  1972  to  July  1974 


Fig  3 Values  of  a(1  — e)  for  Cosmos  462 


Fig  4 Variation  of  perigee  height  y„  for  Cosmos  462 
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Fig  5 Values  of  density  at  a height  of  245  km  in  1972,  P245'  density  index  D , « 
radiation  energy  Fiq  y , Sun-perigee  angle  and  geomagnetic  index  Ap 
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Fig  6 Values  of  density  at  a height  of  240  km  in  1973,  P24o>  ond  density  index  0 , 
radiation  energy  F-|q  7 , Sun-perigee  angle  and  geomagnetic  index  Ap 
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Fig  7 Values  of  density  at  a height  232  km.  P232'  ^ January  to  30  September  1974, 

and  density  index  D , with  solar  radiation  energy  y , Sun-perigee  angle  and 
gaomagnetic  index  Ap  I 
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Fig  8 Values  of  density  at  a height  of  213  km,  P2‘|3>  from  1 October  1974  to  12  March  1975, 

and  density  index  D , with  solar  radiation  energy  F^Q  y , Sun-perigee  angle  and 
geomagnetic  index  Ap 
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nsity  near  the  average  of  both  the  semi-annual  and 
variations,  compared  with  values  from  CIRA  1972 


mean  of  Oct  max /Jan  min  and  Apr  max /Jan  min 


Fig  13  Maximum/minimum  density  ratio  in  the  semi-annual  variation, 
from  1966  to  1975 


